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Abstract 



We show how two important types of phase transition in large N c gauge theory with 
fundamental flavours can be cast into the same classifying framework as the meson- 
melting phase transition. These are quantum fluctuation induced transitions in the pres- 
ence of an external electric field, or a chemical potential for R-charge. The classifying 
framework involves the study of the local geometry of a special D-brane embedding, which 
seeds a self-similar spiral structure in the space of embeddings. The properties of this 
spiral, characterized by a pair of numbers, capture some key universal features of the 
transition. Computing these numbers for these non-thermal cases, we find that these 
transitions are in the same universality class as each other, but have different univer- 
sal features from the thermal case. We present a natural generalization that yields new 
universality classes that may pertain to other types of transition. 



1 Introduction 



A lot of attention has been focused on the properties of the system consisting of the intersection 
of N c color Dp-branes and Nf flavour Dg-branes (p < q). In the large N c ^> Nf limit the Dp- 
branes can be substituted by their corresponding black p-brane supergravity background, while 
the Dg-branes are in the probe limit pp. 

In addition the Dg-branes are extended along q — p of the 9 — p dimensions transverse to the 
Dp-brane, and as a result their gauge degrees of freedom are frozen compared to those of the 
Dp-brane. The dynamics of the p-q strings (which transform in the fundamental of the SU(N C ) 
low energy gauge theory) and the q-q strings are described by the Dirac-Born-Infeld action of 
the Dg-branes. 

Among the issues of interest was the study of the thermodynamic properties of the dual [21 El H] 
Yang-Mills theory and certain thermal phase transitions in the dynamics of the fundamental 
matter. The first study of this nature was for the D3/D7 system and was considered in ref.[5], 
where the authors considered the near-horizon limit of the non-extremal black 3-brane solution, 
corresponding to the AdSs-BHxS 15 geometry (the anti-de Sitter (AdS) spacetime contains a 
black hole). The D7-brane wraps a S 3 C S 5 and extends in the radial direction of the AdSs- 
BH. The size of the S 3 varies as a function of the radial coordinate. The D7-brane embeddings 
then naturally form two classes: embeddings that reach the horizon and hence fall into the black 
hole, and embeddings for which the wrapped S 3 shrinks to zero size at some radial position. For 
these, the D7-brane world-volume simply closes smoothly before the horizon. In an Euclidean 
presentation, the compact, unbounded parts of the D7-brane have the topology S 3 x S 1 since 
the Euclidean time has a periodicity set by the inverse temperature of the system. The classes 
are then distinguished by one or the other compact space shrinking away. The authors of ref. [5] 
proposed that the (topology changing) transition of the D7-brane embeddings corresponds to 
a type of confinement /deconfinement phase transition, now in the meson sector of the theory. 
This system has been extensively studied in refs.[8]-[2"T] and it was shown that it is a first order 
phase transition providing a holographic description of the meson melting phase transition of 
the fundamental matter. 

There is also a unique critical embedding separating those two classes. This solution reaches 
the horizon and has a shrinking S 3 . It has a conical singularity. Solutions of this type will 
occupy much of our attention in this paper. Many of these features generalize to the general 
Dp/Dg system. In ref.[H] the Dp/Dg system was considered and some universal properties, 
associated with this critical solution separating the two classes of embedding, were uncovered. 
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In particular it was shown that for a certain temperature the theory exhibits a discrete self- 
similar behavior, manifested by a double logarithmic spiral in the solution space. This space of 
solutions is parameterized by the bare quark mass and the fermionic condensate. (Geometrically 
these correspond, respectively, to the asymptotic separation of the D7- and D3 branes and 
the degree of bending of the D7-branes away from the D3-branes.) 

The region of solution space where the self-similar spiral is located is unstable, in fact: There 
is a first order phase transition associated with the physics of the system jumping between 
branches of solutions and bypassing it entirely. Nevertheless, it seems that important features 
of the full physical story can be captured by examining the neighbourhood of this critical 
solution. It is remarkable that the critical exponents (or better "scaling exponents" , so as not 
to confuse the physics with the nomenclature of second order phase transitions) characterizing 
this logarithmic structure exhibit universal properties and depend only on the dimension of the 
internal S n wrapped by the Dg-brane. The precise value of the critical temperature is irrelevant. 
The structure is determined by focusing on the local geometry near the conical singularity of the 
critical Dg-brane embedding, and the exponents are then naturally determined by the study 
of possible embeddings in a Rindler space [3 E|. 

The studies described above concern a thermally driven phase transition. As the temperature 
passes a certain threshold, thermal fluctuations seek out the new global minimum that appears 
and the system undergoes a transition to a new phase. In this paper we study transitions of 
the system under the effect of two different types of control parameters: an external electric 
field and an R-charge chemical potential, revisiting work done on these systems in refs. [T6l IT8] . 
We show that the corresponding scaling exponents are again universal and depend only on the 
dimension of the internal sphere wrapped by the Dg-brane. We find that the key properties of 
the critical solution can be determined from the local properties of the geometry, and we find 
that this geometry arises naturally by working in a rotating frame, arrived at using T-duality. 
The resulting physics is not controlled by thermal dynamics, the local geometry is not Rindler, 
and so the exponents are different. The phase transition is driven by the quantum (as opposed 
to thermal) fluctuations of the system, as can be seen from the fact that they persist at zero 
temperature. It is satisfying that we can cast these different types of transition into the same 
classifying framework. 

The structure of the paper is as follows. We begin in section 2 by reviewing the results of 
refs. Hj for the thermally driven phase transition, focusing on the structure of the unstable 
critical solution, extracting the universal properties of the corresponding scaling exponents. We 
highlight the natural appearance of a Rindler geometry. 
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In section 3.1 we consider the case of an external electric field and in the insulator/conductor 
phase transition discussed in ref. [IS]. By employing an appropriate T-dual description of the 
system we demonstrate that the structure of the instability and the scaling exponents can be 
naturally studied by classifying the possible embedding in a flat rotating frame. We observe that 
these scaling exponents are again universal and depend only on the dimension of the internal 
S n sphere, wrapped by the Dg-branes. 

In section 3.2 we consider instead the presence of a finite R-charge chemical potential in the 
Dp/Dq system and demonstrate that the resulting phase transition has the same scaling expo- 
nents as the insulator /conductor phase transition driven by an external electric field. 

We consider some generalizations of the discussion in section 4 and close with some remarks in 
section 5. 



2 Thermal Phase Transition 



Let us begin by reviewing the result of refs.[7J [9]. We will be using the notations of ref. [9]. 
Consider the near-horizon black Dp-brane given by: 



■* ^-fdt 2 + j2dx^j +h* (j^+u 2 dni_ p 



ds 2 = n-- \ - fdt- } j\x: J ■■ /:/- ( — irdn't,,) . <u 

e* = g s H {3 - p)/ * , Obi.., = H- 1 , 

where H(u) = (R/u) 7 ~ p , f(u) = 1 — (uh/u) 7 ~ p and R is a length scale (the AdS radius in the 
p = 3 case). According to the gauge/gravity correspondence, string theory on this background 
is dual to a (p + l)-dimensional gauge theory at finite temperature. Now if we introduce Dg- 
brane probe having d common space-like directions with the Dp-brane, wrapping an internal 
S n C S 8 ~ p and extended along the holographic coordinate u, we will introduce fundamental 
matter to the dual gauge theory that propagates along a (d + l)-dimensional defect. 

If we parameterize S 8 ~ p by: 

dQl_ p = d9 2 + sin 2 Odnl + cos 2 6dn 2 ^ p _ n , (2) 

where dQ^ is the metric on a round unit radius m-sphere, the DBI part of the Lagrangian 
governing the classical embedding of the probe is given byQ 

£ oc e~*J-\g a p\ = —u n sm n 8y/l + fu 2 9' 2 (3) 
9s 



1 We consider only systems T-dual to the D3/D7 one, which imposes the constraint p — d + n + 1 =4. 
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The embeddings split to two classes of different topologies: "Minkowski" embeddings, which 
have a shrinking S n above the vanishing locus (the horizon) and yield the physics of me- 
son states and "black hole" embeddings that reach the vanishing locus, corresponding to a 
melted/deconfined phase of the fundamental matter. These classes are separated by a critical 
embedding with a conical singularity at the vanishing locus, as depicted in figure [T] 




Figure 1: Schematic diagram depicting the Minkowski (left) and black hole (right) embedding so- 
lutions that are separated by a "critical" embedding (centre), which has a conical singularity at the 
event horizon. 

It is convenient to introduce the following coordinates: 

= \ (vT? + \/u 7 -p-u 7 h p ^ , (4) 
L = r cos 9 , and p = r sin 6 . 

Then one can show[T|,[6] that the asymptotic behavior of the embedding at p — ► oo encodes the 
bare quark mass m q = m/27ia' and the quark condensate (4>ip) oc — c of the dual gauge theory 
via the expansion: 

L(p)=m+-^ T + ... (5) 

After solving numerically for each embedding of the Dg-brane, the parameters m and c can 
be read off at infinity. From the full family of embeddings, a plot of the equation of state of 
the system c(m) can be generated. The resulting plot for the D3/D7 system[10] is presented 
in figure [2j The two different colors (and line types) correspond to the two different classes of 
embeddings. The equation of state is a multi-valued function, and there is a first order phase 
transition when the free energies of the uppermost and lowermost branches match. 

The main subject of our discussion is the spiral structure in the solution space near the critical 
embedding [SI E]. In the enlarged portion on the right in figure [2] it is located to the lower left, 
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Figure 2: Plot of the equation of state c(m). The zoomed region shows the location of the first order 
phase transition. There is a spiral structure hidden near the "critical" solution in the neighbourhood 
of m = 0.9185, -c = 0.0225 



roughly at m — 0.9185, — c = 0.0225. The spiral structure that is hidden near this point is a 
signal of the discrete self-similarity of the theory near the critical solution. 

In order to understand the origins of the spiral, we zoom into the space-time region near the 
tip of the cone of the critical embedding[7J E] using the change of variables: 



u = u H + nTz 2 ; 9 = |(^) 



3-P 
4 



f—p 



Here T is the temperature of the background given by: 

7 5 — p 
j -V fu H \ — 

AttR \ RJ 

Leaving only the leading terms in z results in the following metric: 

ds 2 = -(27iT) 2 z 2 dt 2 + dz 2 + dy 2 + y 2 dQ 2 n + dx 2 d + 



(6) 



(7) 



(8) 



The metric (|8j) corresponds to flat space in Rindler coordinates. The embeddings of the Dg- 
branes in the background (j8|) again split into two different classes: Minkowski embeddings 
characterized by shrinking S n (y = 0) at some finite zq, and black hole embeddings, which 
reach the horizon at z = for some finite y = yo (the radius of the induced horizon). The 
equation of motion is derived from the Dirac-Born-Infeld action of the Dg-branes, which has 
the following Lagrangian: 

£ oc zy n y/l+yv . (9) 
The equation of motion derived from this reads: 



zyy" + (yy' -nz)(l + y' 2 ) = . 



(10) 
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Solutions of this equation enjoy the scaling property y(z) — > -y(^z), in the sense that if y(z) 
is a solution to the equation (10) so is -y{^Lz). Under such a re-scaling the initial conditions 
(20,2/0) for the two classes of embeddings scale as: 

z ->z //j,; 2/0 ^W/U; (11) 

This suggests the existence of a critical solution characterized by z = y = 0. One can check 
that y = \fnz is the critical solution. It has a conical singularity at y — z = 0. 



To analyze the parameter space of the solutions we can linearize the equation of motion (10) 
near the critical solution by substituting y(z) = y/nz + £(z), for small £(z). The resulting 
equation of motion is: 

^"W + (« + l)(^)+(W) = 0, (12) 
which has a general solution of the form: 

£(z) = — (A cos(a n In z) + B sin(a n In z)) , (13) 

with r n = ^ ; a n = -y/A(n + 1) - n 2 . 
Note that a n are real only for n < 4, which are the cases naturally realized in string theory 



Now the scaling property of equation (10), combined with the form of the solutions (13) sug- 
gests the following transformation of the parameters (A, B) under the re-scaling of the initial 
conditions given in equation (11): 

(A'\ _ 1 / cos(a n ln/i) sin(a n ln/i)\ / A\ , , 

\B'J ^n+i \^-sin(a„ln//) cos(a n In fx) J \BJ ' ^ ' 

For a fixed choice of the parameters A and B, the parameters (A', B') describe a double spi- 
ral, whose step and periodicity are set by the real and imaginary parts of the critical/scaling 
exponents r n and a n . 



Equation (10) has a Z2 symmetry[7j relating the two classes of solutions (Minkowski and black 
hole embeddings). If the parameters (A,B) describe one class of embeddings, then the pa- 
rameters (—A, —B) describe the other. In this way the full parameter space near the critical 
solution (given by A = 0, B = 0) is a double logarithmic spiral. 

This self-similar structure of the embeddings near the critical solution in our Rindler space is 
transferred by a linear transformation to the structure of the solutions in the (m, c) parameter 
space. If we call (m*, c*) the parameters corresponding to the critical embedding from figure [TJ 
then sufficiently close to the critical embedding we can expand: 

m — m*\ „ , ( A \ , m l2 \ r\i d2 



-Mr B )+0(A') + 0(B') + 0(A,B) . (15) 
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The constant matrix M cannot be determined analytically and depends on the properties of 
the system. Generically it should be invertible (numerically we have verified that it is) and 
therefore in the vicinity of the parameter space close to the critical embedding (m*,c*) there 
is a discrete self-similar structure determined by the transformation: 

m'-m*\_ 1 M f cos(a n ln/i) sin (a n ln/j,)\ (m - m*\ ^ 



d — c* J /i r ' i+1 sin (a n In /i) cos(a n In fj) J \ c — c 

Let us define two solutions to be "similar" if: 

|m' — m"'\ \ 1 / |m — m 



\d -C*| J fjTn+l \ \ c -c*\ J ' ^ 

Then one can see from equation (16) that this is possible only for a discrete set of /is given by: 

fJt=e kn/a n . k = l,2,... (18) 



Note that in general the matrix M in equation ( 16 ) will deform the spiral structure given by 



the transformation (14). However the scaling properties of the theory remain the same as they 
are completely determined by the scaling exponents: r n , a n . Furthermore one can see that 
the scaling exponents depend only on the dimension of the internal sphere S n wrapped by the 
Dg-brane and are thus universal, in the sense that the detailed value of the critical temperature 
is irrelevant. It is the spiral structure that ultimately seeds the multi-valuedness of the space 
of solutions, twisting the (m, — c) curve back on itself as in figure [2j Therefore, it is the spiral 
- and the neighbourhood of the critical solution from where it emanates — that is responsible 
for the presence of a first order phase transition in the system. Whether there is a spiral or 
not can be read off from the scaling parameters (r„, a n ), and since [H] for all consistent Dp/Dq 
systems the condition n < 4 is satisfied the corresponding thermal phase transition (meson 
melting at large iV c ) is a first order one. 



3 Quantum— Induced Phase Transitions 

In this section we will consider a different class of phase transitions. These are arise in the pres- 
ence of external fields, and can happen even at zero temperature, and so since the fluctuations 
driving the transition are no longer thermal, they might be expected to be in a different class. 
Naively, the broad features of the equation of state — multi-valuedness and so forth — have 
similarities with the thermal case, and so it is natural to attempt to trace the extent to which 
these similarities persist. We will find that once we cast these systems in the language of the 
previous section, the similarities and differences will be quite clear. 
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We will first concentrate on the case of an external electric field. The flavoured system, at large 
enough electric field, has an insulator/conductor phase transition, as studied in ref.pi)]- As with 
the thermal transition of the last section, the mesons dissolve into their constituent quarks, but 
this time it is due to the electric field overcoming their binding energy. The transition is of first 
order. 

As we saw in the previous section the scaling properties of the thermally driven phase transi- 
tion are naturally studied in a Rindler frame with a temperature set by the temperature of the 
background. In ref-JTB] it was shown that in analogy to the thermally driven phase transition 
there is a nice geometrical description of the electrically driven phase transition, and the struc- 
ture of the system can be again characterized by an unstable critical embedding with a conical 
singularity at an appropriate vanishing locus (analogous to the event horizon). Here, we will 
generalize this description to the case of the Dp/Dq system. 

Furthermore after an appropriate T-duality transformation we will show that the vanishing 
locus corresponds to an effective "ergosphere" due to a rotation of the coordinate frame along 
the compact directions of the background. The instability near criticality is then naturally 
interpreted as an instability due to the over-spinning of the D(q — 1) brane probes (in the 
T-dual background) as they reach the ergosphere. We then study the structure of the theory 
near criticality by zooming in on the space-time region in the vicinity of the conical singularity. 
Once again, we will find that the structure is entirely controlled by the dimension of the internal 
sphere, S n , wrapped by the D(q — l)-branes (in the T-dual background) — details such as the 
value of the electric field and the temperature of the system, are irrelevant. 



3.1 Criticality and Scaling in an External Electric Field. 

Let us consider the near-horizon black Dp-brane given by the background in equation ([I]). 
Following a similar idea [12] for producing a background magnetic field, if we turn on a pure 
gauge .B-field in the (t, x p ) plane[13j [161 E]> m the dual gauge theory this will correspond to 
an external electric field, oriented along the x p direction: 



B = Edt A dx p . (19) 



The resulting Lagrangian is: 



C oc e -*J-\g aP + B aP \ = -J — ~~~~u n sin" 9y/l + fu*6* . (20) 

9s V / 
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This leads to the existence of a vanishing locus at u = u* given by: 



<~ v = u 7 H p + E 2 L 7 -? , (21) 



at which the action (20) vanishes. Notice that this is distinct from the horizon, and even at 
zero temperature will be present. A study of the local physics near this locus will therefore 
pertain to non-thermal physics. 

The embeddings split into two different classes: Minkowski embeddings which have a shrinking 
S n above the vanishing locus and correspond to meson states and embeddings reaching the 
vanishing locus, corresponding to a deconfined phase of the fundamental matter. These classes 
are separated by a critical embedding with a conical singularity at the vanishing locus. Our goal 
is to explore the self-similar behavior of the theory near this critical embedding and calculate 
the corresponding scaling exponents. 

In order to make the analysis closer to the one performed in refs.[?l |9], for the thermal phase 
transition (described in the last section), we T-dualize along the x p direction. This is equivalent 
to a trading of the pure gauge I?-field for a rotating frame in the T-dual background. Indeed 
the geometry T-dual to equation ([I]), with the 5-field given by equation (19), is given by: 



i ~ i i f d 2 \ 

d~ s 2 = H-^-fdt 2 + dx i) + ZH^Edtdxp + m ( — + u 2 dQ 2 8 __ p + dx 2 p J , (22) 

i=i ^ J ■' 

c* = ^-5; /=!-(- 



i=l 



The background given by equation (22) corresponds to the near-horizon limit of a stack of N c 
D(p — l)-branes smeared along the coordinate x p . Now if we place a probe D(g — l)-brane 
having (d— 1) spatial directions shared with the D(p — 1) — branes, filling the radial direction 
u and wrapping an internal S n inside the S 8 ~ p sphere of the background, we will recover the 



action (20), as we should. 



Note that in these coordinates we have an effective "ergosphere" coinciding with the vanishing 



locus given by equation (21). Now the critical embedding is the one touching the ergosphere 
and having a conical singularity at u = it*. In the (m,c)-plane this embedding corresponds to 
the center of the spiral structure (m*, c*). 

Despite the analogy with the analysis of the thermal phase transition, in this case there is a 
crucial difference, because of the necessity (from charge conservation) for the D(g — l)-brane 
to extend beyond the ergosphere. Indeed since the D(g — l)-brane is an extended object one 
can find static solutions that extend beyond the ergosphere and are non-superluminal. To this 
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end one should allow the D(g — l)-brane to extend along the direction of rotation x p . In the 
original coordinates (before T-dualization) this is equivalent to a non-trivial profile for the A p 
component of the gauge field, which corresponds to the appearance of a global electric current 
along the x p -direction[13, 16J. This is the reason why we refer to the corresponding phase 
transition as an insulator /conductor phase transition. After the transition, the quarks are free 
to flow under the influence of the electric field, forming a current. 

Let us describe how this procedure works in the case of a general D(p — 1)/D(g — l)-intersection. 



Again we will work in the T-dual background (22). Let us consider an ansatz for the D(g — 1)- 
brane embedding of the form: 

6 = 6{u) ; x p = x p (u) ; (23) 

this leads to the action: 




F 2 H I f 2 

' -u n sin n 6Jl + fu 2 8> 2 + —^——x' 2 . (24) 



/ V / — E 2 H 

Now after integrating the equation of motion for x p and plugging the result in the original 
Lagrangian, we get the following on-shell Lagrangian: 




f J ■ 2nn ^ 2n OVT+f^ 2 ■ (25) 

ju Zn sin o — K l 



It is easy to verify that if we choose the integration constant K 2 in equation (25) to satisfy: 

K 2 = E 2 H,u 2n sin 2n # , (26) 



then the action (25) is regular at the ergosphere (u = u*). Note that at the critical embedding 



#o = 0* = and the constant in equation (26) is zero. This constant is proportional to the global 
electric current along the x p direction of the original Dp/Dg-brane system. (See refs. [T3l [T6] 
for a discussion in the case of the D3/D7 system.) 

We are interested in the scaling properties of the theory, near the critical embedding solution. 



Despite the fact that the Lagrangians (20) and (25) describing the Minkowski and ergosphere 
classes of embeddings are different, the fact that at the critical embedding they coincide (K 2 =0) 
shows that the corresponding equations of motion share the same critical solution. Furthermore, 
as we will see, the critical exponents are the same for both types of embedding. 

Let us introduce dimensionless coordinates by the transformation: 
U = U * +Z Y~T V '^ 9= r{ji) * ' Xi (~^) 4 ~* x *> t {~^) 4 ~^ t] H ? E %P^Xp, (27) 
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where D 2 = (7 — p) 2 f*/H*ul, iJ* = H\ u=Ut and /* = f\ u=Ut - To leading order in z and y the 
metric (22) is given by: 

1 P ~ l 

ds 2 = -Dzdt 2 + dz 2 + ci?/ 2 + y 2 dQ 2 n + HiuldQ 2 _ p _ n + 2dtdx p + -^rrdx 2 p + ^ dx 2 . (28) 



i=l 



First consider the case of Minkowski embeddings, characterized by a distance z above the 
ergosphere at which they close (y = y(zo) = 0). The Lagrangian describing the D(g — l)-brane 
embedding is: 

L oc y V/ Vl + y' 2 , (29) 
The corresponding equation of motion is given by: 

Tyl/2 y_ I _ T1 „,n-l_l/2 



9 Z J/ n Z 



a/TT 



,/2 



. 



Equation (30) possesses the scaling symmetry: 



(30) 



(31) 



in the sense that if y = y{z) is a solution to equation (30) so is the function ^y^z). Now under 
the scaling (31) the boundary condition for the Minkowski embedding scales as z — > z /fi. 
This suggests the existence of a limiting critical embedding with z = 0, and indeed: 

y(z) = V2^z , (32) 



is a solution to the equation of motion in equation (30). The corresponding D(g — l)-brane 
has a conical singularity at y = z = 0. Now before we linearize equation (30) and calculate 
the critical exponents let us consider the case of the ergosphere class of solutions characterized 
by the radius of the ergosphere induced on their world-volume. Because of the possibility to 
extend beyond the ergosphere we should consider the analog of the ansatz from equation ([231: 



y = y{z)\ x p = x p (z) 



The corresponding Lagrangian is: 



£* oc y n JDz(l + y' 2 ) + (Fz + 1)£' 2 , 



(33) 



(34) 



where F = D/E 2 H*. After integrating the equation of motion for x p and substituting it into 



the Lagrangian (34), we obtain the following on-shell Lagrangian: 

z^ 2 VFz + ly 2n 



£* oc 



j{Fz + l)y 2n -yl n 



12 



(35) 
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It is easy to see that the Lagrangian (35) is regular at z = 0, y — uq. The equation of motion 



for y(z), derived from the Lagrangian (34) and with the substituted solution for x p (z) is: 



d_ 
dz 



z^y' 



1 + 2/ 



12 



( Fz+1)y 2n_ y 2 

Fz + 1 



2 n. 



ny^z 1 / 2 . 



Fz + 1 



(Fz + l)y 



2n 



v 7 ! + y' 2 = o . (36) 



It is easy to check that equation (32) is a solution to equation p6[). Furthermore for zCl/F 



one can see that equation (36) has the scaling symmetry (31) (note that equation (31 ) suggests 
Vo —* Vol I 1 )- Linearizing equations (30) and (36) near the critical solution (32) by substituting: 



y(z) = V2nz + £(z) 

results in the same equation: 

+ (n + l/2)(z?{z)+t(z)) = 



The general solution of equation (38) is given by: 

£(z) = — (v4cos(o; n ln2;) + B sin (a n In z)) 



where the scaling exponents are given by: 

2n - 1 



a, 



V7 + 20n - An 2 



(37) 



(38) 



(39) 



(40) 



4 7 " 4 

Note that the scaling exponents again, while quite different from those of the thermal case (see 



equation (13)) depend only on the dimension of the internal S n wrapped by the Dg-brane and 
are thus universal for all Dp/Dq systems. Furthermore the discrete self-similarity holds for 
n < 5. By similar reasoning to the thermal case [14], since for all consistent systems realized 
in string theory we have that n < 4, for such systems we may expect that the electrically 
driven confinement /deconfmement phase transition is first order and has the described discrete 
self-similar behavior near the solution that seeds the multi-valuedness of the equation of state. 

The rest of the analysis is completely analogous to the thermal case considered in the previous 
section. Therefore we come to the conclusion that close to the critical embedding (specified 
by m* and c*) the theory has the following scaling property: 



m — m 
c'-c* 



1 M ( cos(a n ln/i) sin (a n In /i)\ j fm-m* 
fx rn+1 \— sin (a n In ji) cos(a n In jj) J y c — c* 



(41) 



with r n and a n given by equation (40) 
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It is interesting to compare the analytic results some numerical studies. Let us consider the 



D3/D7 system. From equation (31) on can see that the variation of the scaling parameter /i in 



equation (41) can be traded for the variation of the boundary conditions of the probe, namely 



zq for Minkowski and yo for ergosphere embeddings. On the other hand, close to the critical 



embedding, the change of coordinates in equation (27) suggests that: 



?o oc yo and uq — oc z 



(42) 



where Uq and 9o are the boundary conditions for the embeddings in the original (not zoomed 
in) background. Note that the parameter uq is related to the constituent quark mass M c 
(in the absence of an electric field) via M c = (it — u h)/ i^na'). 

Close to the critical embedding we have that: 



= («o,in - u*)/(u - «*) and /J = z - m /z 



(43) 



for some fixed boundary conditions w ,i n and 6> 0i i n . Now equation (41 ) suggests that for Minkowski 



embeddings the plot of (m — m*)/(u 



\r n +l 



versus a n ln(w 



should be an harmonic 



function of a ra ln(Mo — m*) with a period 2tt. Similarly for ergosphere embeddings the plot of 
(m — m*)/ 9q" +1 versus a n In Q should be a harmonic function of a n In 9 with a period 2tt. Note 
that the physical meaning of #o can be related to the value of the global electric current (see 



equation (26) and the comment below). 



As can be seen in figure [3| for both types of embeddings the numerical results are in accord 



with equation (41 ) and the analytic results improve deeper into the spiral (large negative values 



on the horizontal axis). Our numerical results confirm that the critical exponents are indeed 
r 3 = 5/4 and = v / 3l/4, as the general analytic results yield. 



m — m„ 

(M0-K,)" 4 



1.5 
1.0 
0.5 

-0.5 
-1.0 
-1.5 



-14 -12 ;-10 -8' 

\ / 



w 





-6 t-4 — j— ln("o-«.) 



■14 /-12 -10\ 



-6 




\m—m 
0.4 
0.2 



,l/0,„ 



-0.2 
-0.-6 



Figure 3: The solid line is a fit with trigonometric functions of period 2tt. The plots confirm that the 
critical exponents of the theory are r 3 = 5/4 and 03 = V31/4. 
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3.2 Criticality and Scaling with R— Charge Chemical Potential. 



Now we study the case when the external parameter is an R-charge chemical potential in the 
dual gauge theory. We will consider the system discussed in ref. [18], where a D7-brane probe 
in the spinning D3-brane geometry [311 [32] was considered. 

The relevant geometry is given by: 



ds 2 = A 1 / 2 [-(HiHMi^fdP + ^ 2 dx 2 + + 

3 

+A -i/2 Hi ^(Rdfa - A\dtf + R 2 dtf) 



(44) 



i=l 

where 



2 4 2 2 3 2 

1 i=i 

with \i\ = sin 9, /i 2 = cos 8 simp, fi 3 = cos 9 cosip- (45) 

Here the parameter uh would be the radius of the event horizon if the angular momentum of 
the geometry was set to zero (g^ = 0). The radius ue of the actual event horizon is determined 
by the largest root of f{u) = 0. The temperature of the background is given by[33j: 

m Ue 2 i 2 i 2 i 2 ( ll ( l2 ( lz\ ^ (2 2\f 2 2\ / ac\ 

= V +qi+q2 + % -^) = 2^u 2 H u E {UE ~ U ^ Ue ~ U2) ' (46) 

where u\ and U2 are the other two roots of f(u) = 0. 



The background (44) has an ergosphere determined by the expression: 

3 



ACH^Hs)- 1 / - J2 ^l(^) 2 = • (47) 



i=i 



Since the background (44) is asymptotically AdSs x S 5 , we can "remove" the ergosphere (47) 



by going to a rotating frame. This is equivalent to gauge shifting A\ from (45) such that 
A't = —fi R + A\. The parameters fi R are set by the condition A' t l \ UE = and hence: 



u 2 H qi 



K u E + q i 

From the behaviour at infinity (u —>■ oo), it is clear that [i l R correspond to the angular velocities 
of the frame along fa. In the dual gauge theory these correspond to having time dependent 
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phases of the adjoint complex scalars or equivalently to R-charge chemical potentials for the 
corresponding scalars [ 



In order to restore some of the symmetry of the metric (44), we will consider the case when 
q 2 = 93. This corresponds to having an S 3 (parameterized by ip, $2, 4>z) inside the deformed S 5 . 
Now if we introduce D7-branes filling the AdS-like part of the geometry and wrapping the S 3 , 
we will add fundamental matter to the gauge theory. Furthermore we are free to rotate the 
D7-branes along <pi and the corresponding angular velocity is interpreted as a time dependent 
phase of the bare quark mas^] (Recall that in introducing D7-branes to the D3-brane system 
we actually add flavours as chiral superfields into the Af = 2 gauge theory). If that phase is the 
same as the phase of the complex adjoint scalar, fi R t, it is equivalent to a R-charge chemical 
potential for both the adjoint scalar and the chiral field. 

On the gravity side of the description this is equivalent to letting the D7-branes have the same 
angular velocity fi R as the rotating frame of the background. Moving to the frame co-rotating 



with the D7-brane corresponds to moving back to the gauge choice for A\ from equation (45). 
The price that we pay is that we again have an ergosphere this time given by: 

A(Wa^)- 1 / - Hi sin 2 9{A]f = 0. (49) 



The possible D7-brane embeddings then naturally split into two classes: Minkowski embeddings 
that have a shrinking S 3 above the ergosphere and ergosphere embeddings which reach the 
ergosphere. These classes are again separated by a critical embedding which has a conical 
singularity at the ergosphere. In analogy to the T-dual description of the previous subsection 
for the external electric field case, the ergosphere embeddings will have to be extended along 
<pi so that they can stay non-space-like beyond the ergospher^J However in this paper we 
are interested in the scaling properties of the theory for parameters (m, c) in the vicinity of 
the critical parameters (m*,c*), corresponding to the critical embedding. As we saw in the 
previous section, modifying the ergosphere class of embeddings so as to be regular at the 
ergosphere does not alter the properties of the theory near the critical solution. In particular 
the scaling exponents characterizing the discrete self-similar behavior of the theory remain the 
same. So henceforth we will focus on the study of the Minkowski type of embeddings. The 
analysis is completely analogous to the one performed in the previous subsection. 

In order to focus on the space-time region close to the conical singularity of the critical em- 

2 We would like to thank A. Karch for pointing this out to us. 

3 Note that in [IB] the ergosphere class of embeddings are not extended along fa. 
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bedding, we consider the change of coordinates: 



where 



UHqi a 71 y (Kt\\ 

M = Merg + -^— COS 6 = - - — , (50) 



M erg = U H- ( ll ( 51 ) 



is the radial coordinate u of tip of the critical embedding or equivalently the 6 = it /2 point of 
the ergosphere. It can be shown that for the values of q2 for which the geometry is not over 
spun (and so has an horizon) the corresponding value of u erg is real. 

After leaving only the leading terms in z and y, we get: 

u 2 

ds 2 /a' = -B x zdt 2 + dz 2 + dy 2 + y 2 dQl - 2q 1 dtd<p 1 + -£dx 2 + R x d§\ , (52) 

R 2 

dti 3 = di; 2 + sin 2 + cos 2 i)d<& ■ D x = ; R\ = ^ R 2 . 



The metric in equation (52) is of the same type as that in equation (28), namely flat space with 



some compact directions in a rotating frame. Therefore the analysis is completely analogous to 



the one for the electric case and hence the scaling exponents are again given by equation (40) 
with n = 3, because the D7-branes are wrapping an internal S 3 : 

r 3 = 5/4 ; a n = V31/A . (53) 

We can again verify this numerically. It is convenient to do this for the single charge case, 
namely q% ^ 0, q2 = q3 = 0. The plot analogous to figure [3] for the electric case, is presented in 
figure |3} The plot represents the variation of the bare quark mass parameter m as a function 
of the initial boundary condition uq — uh, for Minkowski-type embeddings. The parameter 
corresponds again to the bare quark mass for the state corresponding to the critical embedding. 



The good agreement with the result for the critical exponents in equation (40) is clear, and the 
accuracy of the analytic description improves as we go deeper into the spiral (to the left). 

An important observation is that our result does not depend on the values of the R-charges, 
nor the temperature. In fact, this physics persists at zero temperature, such as at extremality 
with all three charges equal q x = q<i = q% = q, or more generally. (Extremality is when Ue — U\ 



or U2, for which T = 0. See equation (46).) The fact that we have the same structure at zero 
temperature (extremal horizon) further confirms that the key properties of the corresponding 
phase transition is indeed driven by the quantum (rather than thermal) fluctuations of the 
system. 
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Figure 4: Plot of the relation between the bare quark mass parameter m and the distance above the 
ergosphere {uq — uh)- The plot is for q = 0.5 in units in which uh = 1. The solid line is a fit with 
trigonometric functions of a period 2ir. The plot confirms that the scaling exponents of the theory are 
r 3 = 5/4 and a 3 = \/31/4. 

4 Criticality and Scaling: Some Generalizations 

In this section we generalize the procedure for the study of the critical behavior employed in 
all three different systems of phase transition (thermal, or in the presence of electric field or 
R-charge chemical potential). This may lay the groundwork for other types of phase transitions 
that may arise in future studies, seeded by spirals with different universal behaviour. 

Note that in all cases there is some vanishing locus. The different classes of Dg-brane's embed- 
dings are being classified with respect to whether they fall into that vanishing locus, or wrap 
an internal S n sphere that is contracting to zero size above the vanishing locus signaling the 
end of the Dg-brane. 

In all cases there is a critical embedding that separates the two classes of embeddings. The 
critical embedding reaches the vanishing locus and has a conical singularity there at some finite 
radius (u* = u H or u erg for the thermal and R-charge cases). 

The main point of the analysis is that after zooming into the space-time region near the conical 
singularity we obtain the metric: 



where D is a non-essential constant. The Dirac-Born-Infeld Lagrangian of the brane is then: 



Note that to extract the key behavior (that we are studying) of this critical embedding (and its 
neighbourhood) there is no need to modify the embeddings which reach the vanishing locus (as 



ds 2 = -Dz k dt 2 + dz 2 + dy 2 + y 2 dn 2 n + . . . , 



(54) 



C oc z k/2 y n y/l +y' 2 . 



(55) 
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we did for the ergosphere class of embeddings) . The critical solution and the linearized equation 
of motion is the same for both classes. Therefore it is sufficient to consider the Minkowski type 



of embeddings and analyze the Lagrangian (|55|). The resulting equation of motion is: 

z k / 2 y n y' 



0, 



12 



(56) 



It is easy to check that equation (56) has the scaling property (31) and the limiting critical 
solution is given by: 

f2n 




Now after the substitution: 



k 



2n 




we obtain the following linearized equation: 

z 2 e(z) + (n + k/2)(zt;'(z)+t;(z)) = 



The general solution of equation ( 59 ) can be written as: 

1 



(Ij {Acos{a ( n k) )lnz) + Bsm{a ( n k) Inz) 

z r„ 



where 



-(*) 



(n + fc/2-l)/2; = \y/l{n + k/2) - (n + k/2 - l) 2 ; 



(57) 
(58) 

(59) 

(60) 
(61) 



are the scaling exponents characterizing the self-similar behavior of the theory. Both being real, 
they control the shape of the spiral which emanates from the critical solution. The oscillatory 
behavior is present for n < 3 + 2y/2 — k/2. For these values of n the theory exhibits a discrete 
self-similarity and the equation of state c = c{m) is a multi-valued function suggesting that 
the corresponding phase transition is a first order one. 

While there is the possibility of complex scaling exponents and hence possibly second order 
phase transitions (if the multi-valuedness goes away when the spiral does), this is not realized 
in the examples that we know from string theory. 

Note that we have k = 2 for a thermal induced phase transition and k = 1 for the quantum 
induced phase transitions that we studied (external electric field and R-charge chemical po- 
tential), arising from the two most natural types of a vanishing locus that one may have: an 
horizon, and an ergosphere. Perhaps other systems will yield different values of k. 
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5 Closing Remarks 



We have succeeded in casting two important types of phase transition (in large N c gauge theory 
with fundamental flavours) into the same classifying framework as the meson-melting phase 
transition. These quantum fluctuation induced transitions (so-called since they persist at zero 
temperature), resulting in the liberation of quarks from being bound into mesons as a result 
of the application of an external electric field, or a chemical potential for R-charge, turn out 
to have the same underlying structure. It is distinct from that found for thermal fluctuation 
induced transitions. The structures are controlled by the local geometry of the spacetime seen 
by a critical D-brane embedding (it is the borderline case between two physically distinct classes 
of embedding), and while it is Rindler for the thermal case with an horizon at the origin, it is 
(after a T-duality in order to geometrize the discussion as much as possible) a rotating space 
with a simple "ergosphere" type locus. The technique of characterizing the physics in terms 
of this underlying classifying space [TJ is rather pleasing in its utility, and we extended our 
analysis to the natural generalization of this space, extracting the scaling exponents that might 
pertain to physics from future studies. 

Of course, there is much interest in how much we can learn about finite N c physics (for ap- 
plications to systems such as QCD) by studying universal features of large N c . Unfortunately, 
it is almost certain that much of this is far from robust against 1/N C corrections. The spiral 
structure is rather delicate, and the stringy corrections arising in going away from the large N c 
limit would generically severely modify the classifying spacetimes we've been studying, erasing 
the spiral and its self-similarity. The absence of the spiral is necessary for there to be (at best) 
a second-order transition at finite N c , since it results in multi-valuedness of the solution space, 
requiring the system to perform a first order jump. 

It is tempting to speculate, however, that the nature by which the spiral is destroyed by l/N c 
corrections might (especially since the setting is so geometrical) be characterizable in a way 
that allows universal properties of the second (or higher) order phase transitions to be deduced 
from the properties of the spiral at large N c . We leave such explorations for later work. 
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